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Abstract 

In this paper we consider a model of an atom with n energy levels interacting 
with n(n-l)/2 external (laser) fields which is a natm'al extension of two level system, 
and assume the rotating wave approximation (RWA) from the beginning. We revisit 
some construction of analytical solutions (which correspond to Rabi oscillations) of 
the model in the general case and examine it in detail in the case of three level 
system. 
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1 Introduction 



Our purpose of this paper is to apply a method developed in this one to Quantum Com- 
putation. Quantum Computation (or Computer) is a challenging task in this century for 
not only physicists but also mathematicians. See for example jT] as a general introduction. 

Quantum Computation is in a usual understanding based on qubits which are based 
on two level system (two energy levels or fundamental spins of atoms), See j2j, 
jS] as for general theory of two level system. In a realistic image of Quantum Computer 
we need at least one hundred atoms. However then we meet a very severe problem called 
Decoherence which will destroy a superposition of quantum states in the process of unitary 
evolution of our system. See for example 0| or recent JU] as an introduction. At the 
present it is not easy to control Decoherence. 

By the way, an atom has in general infinitely many energy levels, while in a qubit 
method we use only two energy levels. We should use this possibility to reduce a number 
of atoms. We use n energy levels from the ground state (it is not realistic to take all 
energy levels into consideration at the same time). We call this n level system a qudit 
theory, see for example [T3j, [THj, [in]- 

In a qubit case we need Rabi oscillations to construct quantum logic gates. See [12] 
as a simple introduction to quantum logic gates. Similarly we also need Rabi oscillations 
to construct them in a qudit space. However a general theory of Rabi oscillations in a 
level system has not been developed enough as far as we know. For three level system see 
[H], [5 a-iid [S]. Therefore we attempt such a theory in this paper. 

In this model we assume the rotating wave approximation (RWA) from the beginning, 
otherwise we cannot solve the model. However there is no problem on the approximation 
in the weak coupling regime (the coupling constant is small enough compared to the energy 
differences of an atom) . We note that there is a problem to use this approximation in the 
strong coupling regime. 

Before stating our result let us make a review of Rabi oscillations (or coherent os- 
cillations) in a two level system. For a quantum version of the Rabi oscillation see for 
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example, 0, [ISl- 

Let {(Ji, (72, (T3} be Pauli matrices : 



1 \ -i \ / 1 

^1=1 ' ^2 = ' ^3 = I , (1) 

10/ \ i \ -1 



and we set 



1 I 1 1 1 

^ I ' ^ 







V 



1 

Let us consider an atom with 2 energy levels E'o and Ei { Ei > Eq). Its Hamiltonian 
is in the diagonal form given as 

' Eq 

^^0=1 |. (2) 



El 



This is rewritten as 



, 1 \ / \ A 

Ho = Eo \ ^ I + (El - Eo) I I = E0I2 + y (I2 - t^s) , 

where A = i?i — E'o is an energy difference. Since we usually take no interest in constant 
terms, we can set 

H, = ^ (I2 - a,) . (3) 

We consider an atom with two energy levels which interacts with external (periodic) 
field with gcos{ujt + 0). In the following we set h = 1 for simplicity. The Hamiltonian in 
the dipole approximation is given by 

H = Hq + gcos{ujt + (f))ai = y (I2 - 0-3) + gcos{ujt + 0)(Ti, (4) 

where u is the frequency of the external field, g the coupling constant between the external 
field and the atom. This model is comphcated enough, see [Zj, [TT], (THj, [IZj, [TH] . 

In the following we set = for simplicity and assume the rotating wave approxima- 
tion (which neglects the fast oscillating terms), namely 

cos(cut) = ^(e'"* + e"'"*) = ^e''"*(l + e""'"*) ^ Je'"*, 

Zi 



and 



cos(a;t)ai 



cos{ujt) 



y cos{ujt) 
therefore the Hamiltonian is given by 



1 
2 



e 



ILdt 



^ = f (l2 - ^3) + f (e^-* 



e (T+ + e a_ 



A , 

Y (I2 - CT3 



) + i7 (e^"V+ 



by the redefinition of (7 {g/2 — > g). It is exphcitly 



H 



ge 



We would hke to solve the Schrodinger equation 



dt 



For that let us decompose H in ^ into 



^e*^* 



1 



^? 
^7 A 



SO if we set 



-Auit 



-iujt 



then it is not difficult to see 



(it 



/ 



^7 



V ^ A- w 

which is easily solved. For simplicity we set the resonance condition 



A = a;, 



then the solution of (fTIH) is 



1 



^t) = exp { -igt I I ) $(0) 

1 



cos{gt) —ism{gt) 
-ism(gt) cos(gt) 



As a result, the solution of the equation ((Tj) is given as 



I 



\ 



-iujt 



( 



\ 



I 



-iuit 



cos{gt) —ism{gt) 



^ —ism{gt) cos{gt) 



m (12) 



by ©• If we choose $(0) = (1, 0)^ as an initial condition, then 



/ 



cos{gt) 



y —ie ^'^*'sm{gt) 

This is a well-known model of the Rabi oscillation (or coherent oscillation). 
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2 N Level System • • • General Theory 

In general, an atom has an infinitely many energy levels, however it is not realistic to 
consider all of them at the same time. Therefore we take only n energy levels from the 
ground state into consideration (-En-i > ■ ■ ■ > -Ei > Eq). Then from the lesson in the 
preceding section the energy Hamiltonian can be written as 







Ai 



(14) 



where = Ej — Eq for 1 < j < n — 1. For this atom we consider the interaction with 
n{n — l)/2 independent external (laser) fields corresponding to every energy difference 
{Ej — Ei for j > i). Then the interaction term assuming the RWA from the beginning is 
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given as 



V 



-iwQ2t 



-iu)i2t 



.n-lt \ 



-^iuJ2.n-lt 



-iWn-2,n-lt g 



SO the full Hamiltonian is 



H = Ho + gV, 



where g is the coupling constant. Here we set 



^1 — 1^01) ^2 — ^12, ■ ■ ■ ; ^n-2 — ^n-3,n-2, — ^n-2,n-l 



for simplicity. We would like to solve the Schrodinger equation 



dt 



Similarly in (jH} let us decompose H : For 



U 



J,{u)-i+W2)t 



-^i{tjJ-i+uJ2^ hi^n-2)i 



-^i{uJi+U)2-\ \-<^n-l)t 
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it is easy to see 



Hu = U^HU 



where 

for j — i > 2. 
By setting 








g 


gQieo2t 


gQi^0,n-2t 


^g^^O,n— 1 1 
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Ai 


g 


gQiei,n-2t 




gQ — i€02t 




g 




gQie2,n-2t 




gQ—i€0,n-2t 




-iei,n-2t 


gQ—i€2,n-2t 


An-2 


g 


g^ i^O,n—li 






gQ—i€2,n-lt 


g 


A„-i 








+ OUi+2 H \-UJj) 



















(19) 



it is not difficult to see 



where 



at 



(20) 



Hu 
( 



g^ 





g 

—ieo2t 



Ai - a;i g 

g A2 - {uji + UJ2) 



gQ—ieo,n-2t gQ—iei,n-2t 

gQ ^^0,71 — 1^ g^ 1'^l,n — li 



ge 



ge 



i^0,n-2t 



ge 



«ei,n-2t 



ge 



i^2,n-2t 



ge-ie2,n-2t . . An-2 - ^2?=! 

-«e2,n-lt 



ge 



«ei,n-it 



^e 



ii2,n-lt 



g 



g 



An_i — Xir=l '^l J 



(21) 
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At this stage we set the resonance conditions 

Ai = a;i, A2 = a;i + a;2, • • • , A„_2 = ^ uji, A^-i = J2^i 



n-2 



n— 1 



1=1 



1=1 



Uj = Ej - Ej_i (j = 1, 2, • • • , n-1) 



to make the situation simpler. Therefore Hu can be written as 

Hu = 9C + gR, 



where 



C 



1 

1 1 

1 1 



1 1 
1 



and 













g*eO,n-lt 











g«ei,n-2t 


giei,n-lt 


Q—ieo2t 








Qie2,n-2t 


gie2,n-lt 



v 



Q—ie0,n-2t g— jei,n-2t g— 162, 71-2* 

g — iCQ.n-lt g — *<:i,n-lt Q — i^2,n-lt 

Therefore the Schrodinger equation becomes 

zj^^ = {gC + gR)^. 













Now it is easy to solve 



which solution is 



dt 



^{t) = exp{-igtC)^{0). 
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(22) 



(23) 



(24) 



(25) 



By making use of the method of constant variation 



^(t) = exp(-i^tC)^(0) — > = exp{-igtC)(p{t) 



and substituting this into we have the equation 
.d 



dt 



gexp{igtC)R{t)exp{—igtC)(j) = gA{t)(f), where 



A{t) = exp{igtC)R{t)exp{-igtC). 
For this equation a formal solution called the Dyson series^ has been known 
(f){t) = : exp (^-ig A{s)ds^ : (f)(0) 



(26) 
(27) 



1-igJ A{s)ds - g"^ J |v4(s) J A{u)duj ds + 



0(0). 



(2^ 



However it is usually almost impossible to perform integrals. Anyway, our formal solution 
of the equation (pUj) becomes 



^(t) = exp{-igtC) : exp i^-ig A{s)dsj : ^ (29) 
where \1/ is a constant vector. As a conclusion, a formal solution of (|T7jl is given by 

^(t) = U{ty-^{t) = U{tyexp{-igtC) : exp (^-ig A{s)ds^ : * (30) 
with U{t) in (|TR|l . 

Lastly let us calculate exp{—igtC) in ()30|1 . For that purpose we need to make C in 
a diagonal form. First let us calculate eigenvalues. The characteristic equation of C 



A -1 
-1 A -1 

-1 A -1 



IS 



/.(A) = 



-1 A -1 
-1 A 



(31) 



^In Mathematics it was called the Peano series a long time ago, [20! 



By the Laplace expansion with respect to the first column we have a recurrent equation 

/„(A) = A/„_i(A) - /„_2(A) and /2(A) = A^ - 1, /3(A) = A=^ - 2A. (32) 

For example, /4(A) = A"^ — 3A^ + 1 . From these equations it is not easy to conjecture the 
eigenvalues. However we know a good method ! We set 

X^x+-^x + x-\ (33) 

X 

Then for example, 

f2{x) = + 1 + x'"^, f3{x) =x^ + x'^ + x~^ + x"^, /4(a;) = + + 1 + a;"^ + x"^. 
In general 

Ux) = = x-n ^ ^_ (34) 

The solutions of fn{x) = are {exp(v^^7rj'/ (n + 1)) | j = 1, 2, ■ ■ ■ , n, n + 2, ■ ■ ■ , 2n + 1} 
— i) because x 7^ ±1. Therefore our solutions are 

\j — Xj + xj^ = 2cos ^ — I^) J = 1) 2, • • • , n (35) 

by removing the overlap. 

Next let us calculate an eigenfunction corresponding to Xj 

C|j) = A,|j) = 2cos(^) \j) 

for j = 1, 2, ■ • ■ , n. If we define 



T 



' sin f ^!^y sin f sin ( ),■■■, sin f^^) 1 , (36) 



yn + ly \n + lj \n + lj \n + lj \n+l 

then it is not difficult to see 

CU) = 2cos |j) and m ^ 5j,. (37) 

We leave the proof to the readers. Now we define 



O^iOjk) where Ojk ^ ^ ^^^in (^-^j =^ OO^ ^ In (38) 
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and 

D = (Djk) where Djk = 2cos ( ^"^ ) 6jk, (39) 
then we can diagonahze C like 

C = ODO^. (40) 

Therefore we finally obtain 

exp{-igtC) = Oexp{-igtD)0^ (41) 

or in components 

exp^gtC),, = E -P ) sin sin . (42) 

In the case of n = 3 let us write down exp{—igtC) explicitly for the sake of latter 
convenience. 

^ 1 + cos(y2^t) -iV2sm{V2gt) -1 + cos{V2gt) ^ 

-iV2sm{V2gt) 2cos{V2gt) -iV2sm{V2gt) ■ (43) 
\ -1 + cos{V2gt) -iV2sm{V2gt) 1 + cos{V2gt) y 



exp{-igtC) = i 



For (PT|) with (jl^ we would like to calculate the Dyson series fl28|) . However it is 
almost impossible in the general case, so in the last section we calculate the next-leading 
term in the case of n = 3 by making use of 

3 N Level System • • • An Exact Solution 

In this section we present an exact solution in the general case under some "consistency 
condition" (in our terminology). 

In section 2 we consider a very special case : namely, in we take 



eij = ujij - (tUj+i + u;i+2 H \- ujj) = (44) 
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for all j — i > 2. We call this a consistency condition. Then Hu becomes a constant 
matrix ! 

1 1 ■ ■ 1 1 ^ 

1 1 ■ ■ 1 1 
1 1 ■ ■ 1 1 



Hu = gC + gR = g 



1 1 1 ■ ■ 1 
1 1 1 ■ ■ 1 



so it is easy to solve the equation (^3)1 



d ~ 



The solution is formally 



^{t) = exp{-igtQ)^{0). 
Let us calculate exp{—igtQ) explicitly. If we define 

|l) = (l,l,---,l,lf, 

then it is easy to see Q = |1)(1| — !«, so 

^(t) = e'3^exp{-igt\l){l\)^{0). 



(45) 



Since (1|1) = n, 



so that 



{\l){l\f = n'-'\l){l\, 



[-igtf 



igtf 



exp(-z^t|l)(l|) = 1„ + E = In + E 



fc=i 



k=l 



k\ 



n [^^Q fc! J n 



(46) 
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Therefore we obtain 



i(t) = ,.».|i„ + !^EtHwO^|i)(i|U(o). 



n 



(47) 



Here if we choose '^(0) = (1, 0, ■ ■ ■ , 0, 0)^ = |0), then 



n 



As a result, the solution that is looking for is just 



n 



(4J 



with |T8|) . In particular, in the case of n = 3 the solution is 



V 



\ / exp(-i3gf)+2 \ 
3 



-i(ijJi+LLi2)t 



exp(-r3gf)-l 
3 

exp(^i3gf)-l 
3 



igtexp{-ag£)+2 \ 

^ 3 
t(o-^,)f exp(-i3gf)-l 



^i(q-t^-i-t^^)f exp(-i3gf)-l 
3 



(49) 



4 Three Level System • • • An Approximate Solution 

We calculate the case of n = 3 in more detail. For that purpose let us list f/(t)^ 
exp(— i(yftC) and R{t) once more : 

f/(t)t = 



-iuilt 



V 



-i(a;i+LJ2)t 



and 



exp( — ZfyftC) 



1 + cos(V2^t) -iV2sin(V2^t) -1 + cos(V2^t) 



-1 + cos{y/2gt) -iy/2sm{y/2gt) 1 + cos{y/2gt) J 
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and 



exp(iet) ^ 



R{t) = 

y exp(— iet) y 
where e = eo2 = ^0102 — {uji + UJ2) for simplicity. 

First of all we must calculate A{t) in (jSZj)- This calculation is tedious even in the case 
of n = 3. The result is 



A{t) = exp{igtC)R{t)exp{-igtC) 



^ ail ^12 Ois ^ 



0-21 0,22 ^23 
\ C^Sl '^32 '^33 J 



(50) 



where 



ail = — sin^(-\/25ft)cos(et), 

ai2 = v^sin(v^5't)sin(et) — i\/2sm[\/2gt)cos{-\/2gt)cos{et) ^ 
ai3 = {1 + cos^(v^5't)^ cos(et) + 22cos(A/2(7t)sin(et), 

021 = "\/2sin(-\/2(7t)sin(et) + iA/2sin(-\/25ft)cos(-\/2(7t)cos(et), 

022 = 2sin^(v^5't)cos(et), 

a23 = — A/2sin(-\/25ft)sin(et) + i\/2sm{\/2gt)cos{\/2gt)cos{et) , 
(231 = (1 + cos^{\/2gt)^ cos(et) — ■i2cos(-\/2(7t)sin(et), 
(232 = — v^sin(v^5f)f:)sin(e)f:) — 'iv^sin(A/25ft)cos(A/25't)cos(eit), 
%3 = — sin^(-\/25ft)cos(et). 



Next we must calculate the Dyson series ()28j) with A{t) above. To calculate all terms 

in (pH|l is of course impossible, so we calculate up to the next-leading term : 

rt 
10 

with 0(0) being the ground state 0(0) = (1, 0, 0)"^. By the way, it is convenient for us to 
calculate 



0(t) 



1-ig [ A{s)ds 
Jo 



0(0) 



(51) 



exp{—igtC)(f>(t) 



exp{—igtC) — ig exp{ig{s — t)C)R{s)exp{—igsC)ds 
Jo 



0(0) (52) 
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rather than ()5H) itself. The result is 

exp{-igtC)(l){t) = X2(t), X3(t))^ , 

where 



xiit) 



^ + - |(-2 + cos{V2gt)) ^ - 

ig j sm{^/2gt) + sm{{V2g + e)t) sm{V2gt) + sm{{y^g - e)t) 
"8 \ 2V2g + e ^ 2V2g - e 



+ 



X2{t) 



V2g + 



V2g- 



V2 4 
i\/2fl' J sin(et) + sm(\/2gt) sin(et) — sm(^/2gt) 1 

V2g + e ^ V2g-e J ~ 
y/2g j cos{{^/2g + e)t) - cos{V2gt) cos{{V2g - e)t) - cos{V2gt) 1 
~8~ \ 2y/2g + e ^ 2^2^ - e J 



2V2g + e 2^g - e 

-l + cos(v/2,t) _ I ^2 + cos(y2,t)) ^ - 

( sin(v/2^t) + sin((^^ + e)t) ^ sin(v/2^t) + sin((x/2^ - e)t) | ^ 

-cos(^) + cos((^ + e)t) -cos( j|) + cos((^ - e)t) 

V2g-t 



9 I gt . 

2 



(53) 



V2g + t 

We leave the check of this calculation to the readers. 

As a result we obtain the approximate solution (up to the next-leading) in the case 
of n = 3 

^ xiit) ^ 

-'^^'x2it) (54) 









^ Xl{t) ^ 










X2(t) 






^ g-i(iJi+a)2)t ^ 




^ X3{t) ) 





V 



xsit) J 



with xi{t) , X2{t) , X3{t) above. This is our main result. 
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5 Discussion 



In this paper we developed a general theory of the n level system by assuming the rotating 
wave approximation and gave an analytical solution under the "consistency condition". 
Moreover in the three level system we gave some approximate solution (up to the next- 
leading term), which may be enough because we are treating the weak coupling regime 
[g is small enough compared to the energy differences of the atom). To continue further 
calculations (higher order corrections) is not easy at the present. 

Here we state our motivation once more. In Quantum Computation with three level 
system we must first of all construct or approximate all elements in U{3) by making use 
of (a kind of) Rabi oscillations ■ ■ ■ a universality of one qudit theory. This is not an easy 
problem. In fact, in almost all papers in qudit theory this has not been mentioned as 
far as we know. Our paper is an attempt to this problem and we gave the approximate 
solutions (jHH). We need a furthermore study. 

We list the equation in the three level system once more 



.d_ 

'dt 



^2 







-ILOlt 



Ai 



ge 

A2 



^1 



(55) 



and present 



Problem Show a good method to solve (or to solve approximately) the equation above. 

This problem seems to be very hard. 
Moreover, study the general equation in the three level system (without the rotating wave 
approximation) 



.d 

I— 
dt 



^2 

This is our next target 







gcos{iuJit) gcos{iuJo2t) 



\ 



gcos{~iLJit) Ai gcos{iuj2t) 

y gcos{—iu!o2t) gcos{—iu!2t) A2 



^2 



(56) 
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